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FOREWORD 


The  theory  of  the  electronic  structure  of  molecular 
systems  is  discussed.  The  potential  energy  surface  is  described 
as  a  conceptual  and  computational  tool  in  the  study  of  the 
kinetics  of  fast  reactions  which  might  occur  in  detonations. 

It  is  suggested  that  certain  excited  states  are  diminished 
in  energy  during  compression,  and  that  intersections  of 
ground  and  excited  -state  potential  energy  surfaces  occur, 
thereby  creating  much  lower  activation  energy  barriers  and 
greatly  enhancing  reaction  rates. 

Ab  ini tio  and  semiempi rical  calculations  of  ground 
and  excited  electronic  states  are  performed  on  the  prototype 
molecular  explosive  system,  (N0+)n.  Each  N0+  molecule  is 
allowed  to  interact  with  the  other  molecules,  as  a  model 
of  various  states  of  compression.  The  results  of  the  calculations 
show  that  the  intersections  are  forced  to  occur  under  compression, 
and  that  the  activation  energy  for  the  N-0  bond  scission 
process  is  substantially  reduced  from  its  gaseous  state 
value  of  254  kcal/mole.  The  results  also  show  that  IR  and  UV 
spectra  are  shifted  to  lower  frequencies.  These  calculations 
are  the  first  in  a  series  to  be  performed  on  explosive  molecules, 
from  which  criteria  for  sensitivity  to  shock  initiation  leading 
to  detonation  will  be  deduced.  The  author  thanks  Drs.  D.  J.  Pastine 
M.  J.  Kamlet,  E.  T.  Toton,  and  S.  J.  Jacobs  for  helpful  discussions. 
This  work  was  carried  out  under  Task  Number  ZR01301,  IR-221,  and 
Task  Number  ZR01305,  IR-214. 


1/2 


■ 


dtt 


NSWC  TR  79-175 


CONTENTS 


Page 

INTRODUCTION  .  4 

THE  MICROSCOPIC  AND  MACROSCOPIC  ASPECTS  OF  SENSITIVITY 

TO  DETONATION .  6 

THE  ACTIVATION  ENERGY  AS  A  POTENTIAL  ENERGY  CONCEPT  ...  8 

THE  QUANTUM  MECHANICAL  PROBLEM  .  10 

CALCULATIONS  OF  GROUND  AND  EXCITED  STATES  FOR  THE  (NO+)n 

MOLECULAR  SYSTEM  .  26 

SUMMARY  AND  FUTURE  WORK  TO  BE  PERFORMED .  41 


3 


NSWC  TR  79-175 


ILLUSTRATIONS 

Figure  Page 

1  Avoided  Crossing  and  Intersection  for  Two 

Potential  Energy  Curves  .  22 

2  Excitation  Energies  for  (N0+)n  .  28 

3  Excitation  Energies  per  Molecule  for  (N0+)n .  30 

4  Energy  Curves  for  NO+ . 32 

5  Two  Interacting  N0+  Molecules . 33 

6  Energy  Curves  for  Broadside  Collision  .  34 

7  Energy  Curves  for  Head-on  Collision  .  35 

8  Energy  Curves  Resulting  from  Coupling  of 

Electronic  and  Nuclear  Motions . 37 

9  Repulsion  Energies  and  Pg  Values  for  (NO+)2  ...  40 


4 


NSWC  TR  79-175 


INTRODUCTION 


The  sensitivity  of  an  explosive  is  often  categorized  in  terms 
of  its  response  to  stimuli  such  as  impact, ^  friction,  and  shock. ’ 

In  Reference  1,  it  was  demonstrated  for  a  variety  of  high  (secondary) 
explosives  that  an  impact  explosion  behaves  like  a  low  temperature 
thermal  decomposition  which  implies  that  temperature  is  required 
for  initiation  here.  In  that  work,  impact  sensitivity  trends  were 
deduced  on  the  basis  of  molecular  structure.  In  the  case  of  friction 
sensitivity,  temperature  again  determines  initiation  to  explosion. 

As  indicated  in  Reference  3,  the  critical  pressure-time  profiles 

required  to  achieve  a  shock  transition  to  detonation 
are  determined  by  the  chemical  and  physical  nature  of  the  explosive 
as  well  as  the  explosive  dimensions.  Our  knowledge  of  the  chemical 
and  physical  nature  of  an  explosive  derives  from  theoretical  and 
empirical  considerations.  In  view  of  the  complexity  of  the  response 
of  explosives  to  the  various  stimuli,  any  approach  which  attempts 
to  relate  the  totality  of  the  chemical  and  physical  characteristics 
of  an  explosive  to  sensitivity  must  be  semiempir ical  in  nature. 

Pressure  plays  a  crucial  role  in  the  shock  initiation  of 
detonation.  The  pressure-time  profile  relationship  Pnt  =  const, 
has  been  deduced  from  a  theoretical  consideration  of  relaxation 
processes  which  channel  energy  into  an  explosive  molecule. ^  Here, 
the  pressure  P  consists  of  a  temperature-dependent  term  PT  as  well 
as  a  term  arising  from  molecular  interactions  at  0°K,  PQ  =  -dE/dV, 

where  E  is  the  internal  energy  and  V  is  the  volume.  These  contribu¬ 
tions  to  P  are  well  understood  when  examined  from  the  point  of  view 
of  an  equation  of  state.  On  the  other  hand,  the  relative  importance 
of  these  terms  is  not  understood  with  respect  to  the  shock  initiation 
of  detonation  problem.  It  is  well  known  that  shock  compression 
produces  high  temperatures  which  always  accelerate  chemical  reactions, 
and  which  generate  the  term  PT.  For  example,  a  shock  wave  in  liquid 


iRamlet,  M.  J.  ,  ''The  Relationship  of  Impact  Sensitivity  with 
Structure  of  Organic  High  Explosives,  I.  Polynitroaliphatic  Explosives" 
in  Proceedings  of  the  Sixth  Symposium  (International)  on  Detonation, 

24  August  1976,  pp,  312-322, 

2Bowden,  F.  P.,  and  Yoffe,  A.  D. ,  Initiation  and  Growth  of 
Explosions  in  Liquids  and  Solids,  (Cambridge  University  Press,  1952) , 

^Price,  D. ,  "Shock  Sensitivity,  A  Property  of  Many  Aspects,"  in 
Proceedings  of  the  Fifth  Symposium  (International)  on  Detonation, 

1$  August  1970,  pp.  207-217, 

^Pastine,  D.  J,,  Edwards,  D.  J. ,  Jones,  H.  D.,  Richmond,  C.  T. 
and  Kim,  K,,  "Some  New  Concepts  Relating  to  the  Initiation  and  Failure 
of  Detonable  Explosives,"  in  Timmerhaus,  K.  D.,  and  Barber,  M.  S., 
eds.,  High-Pressure  Science  and  technology,  Vol.  2,  (New  Yorks 
Plenum  Publishing  Corp.  1979) ,  pp.  364-382. 

5pastine,  D.  J.,  and  Bernecker,  R.  R. ,  "P,V,E,T  Equation  of 
State  for  1, 3, 5-triamino-2, 4, 6-trinitrobenzene,"  Journal  of  Applied 
Physics,  Vol.  45,  No.  10,  1974,  p.  4458. 
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nitromethane  produces  detonation  within  a  few  microseconds  when  the 
explosive  is  compressed  to  60%  of  its  original  volume  and  the 
temperature  reaches  1140°K.  However,  little  is  currently  known 
about  the  detailed  mechanism  which  relates  Pg  to  initiation  and 
to  the  steady-state  propagation.  This  information  is  needed  for 
a  fundamental  understanding  of  sensitivity  to  detonation  as  well 
as  for  the  correct  classification  of  explosives  according  to  their 
sensitivity,  which  is  essential  for  their  safety  and  handling. 

The  interpretation  of  the  mechanism  of  initiation  requires 
the  consideration  of  all  of  the  microscopic  and  macroscopic 
characteristics  of  an  explosive,  the  former  being  the  most  difficult, 
if  not  impossible,  to  analyze  experimentally  at  high  pressures  and 
temperatures.'  An  understanding  of  the  microscopic  aspects  is  of 
prime  importance  to  the  chemist  who  is  confronted  with  synthesizing 
explosive  molecules.  Conditions  in  the  laboratory  permit  the 
tailoring  of  molecules  which  meet  specific  requirements  for  weapons 
systems.  If  new  insensitive  high  explosive  molecules  are  to  be 
synthesized,  it  should  be  possible  to  recognize  those  pertinent 
intrinsic  structural  features  which  define  sensitivity,  as  distinct 
from  those  apparent  features  derived  from  macroscopic  consider¬ 
ations.  In  this  report,  a  program  is  introduced,  the  ultimate 
objective  of  which  is  to  define  criteria  for  initiation  of  detonation 
and  to  interpret  the  steady-state  detonation  based  on  the  electronic 
structure  of  molecules  comprising  explosives. 

Semiempirical  computational  methods  have  been  used  to  relate 
electronic  structure  and  sensit ivity.8  These  very  approximate 
calculations  dealt  with  only  the  ground  electronic  states.  However, 
a  recent  series  of  calculations  has  been  used  to  determine  a  striking 
correlation  between  sensitivities  to  detonation  and  electronic 
charge  distributions  in  the  excited  states  of  the  molecules  TATB 
and  five  other  nitroaromatics ,  tetryl  and  six  other  nitroamines, 
and  nitroglycerine  and  seven  other  nitric  esters.9  These  calculations 
were  performed  on  the  isolated  molecules,  and  did  not  address  the 
question  of  whether  the  energy  in  the  explosive  is  sufficiently 


^Campbell,  A.  W. ,  Davis,  W.  C. ,  and  Travis,  J.  R. ,  "Shock 
Initiation  of  Detonation  in  Liquid  Explosives,"  Physics  of  Fluids, 

Vol .  4,  No.  4,  1961,  p.  498. 

7Bardo,  R.  D. ,  Hall,  T.  N.,  and  Kamlet,  M.  J. ,  "The  Volume  of 
Activation  in  the  Shock  Initiation  of  Explosives,"  Combustion  and 
Flame,  Vol.  35,  1979,  p.  259. 

“Haskins,  P.  J. ,  "Electronic  Structure  of  Some  Explosives  and 
Its  Relationship  to  Sensitivity,"  Research  in  Primary  Explosives-ERDE 
1975,  Presentation  No.  6,  1975. 

9Delpuech,  A.,  and  Cherville,  J.,  "Relation  entre  la  Structure 
et  la  Sensibility  au  Choc  des  Explosifs  Secondaires  Nitrys-Crityre 
Moiyculaire  de  Sensibility.  I.  Cas  des  Nitroaromatiques  et  des 
Nitramines,"  Propellants  and  Explosives,  Vol.  3,  1978,  p.  169. 
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large  and  localized  to  reach  excited  states.  In  the  present  report, 
we  are  interested  in  interacting  molecules  as  a  simulation  of 
compressed  states,  in  order  to  determine  trends  and  changes  in  the 
ground  and  excited  electronic  states  and  the  corresponding  changes 
in  electronic  excitation  energies.  The  calculations  will  require 
the  development  of  new  theory  based  on  the  wave  equations  of 
Reference  10  and  a  utilization  of  both  ab  initio  and  semiempirical 
computational  methods. 

It  will  be  shown  from  the  methods  of  quantum  mechanics  how 
activation  energies  are  related  to  electronic  structure,  and  how 
these  energies  and  the  corresponding  reaction  rates  change  as  the 
molecular  electronic  structure  is  altered  under  compression.  It 
will  be  seen  that  direct  excitation  of  electronic  states  need  not 
occur  for  certain  excited  states  to  contribute  to  reaction  processes. 
In  the  case  of  interacting  NO+  molecules,  a  0°K  isotherm  will  be 
presented.  The  results  of  preliminary  calculations  on  the  NO+ 
molecules  and  the  relationship  to  sensitivity  will  be  discussed. 


THE  MICROSCOPIC  AND  MACROSCOPIC  ASPECTS  OF  SENSITIVITY  TO  DETONATION 


A  complete  analysis  of  the  sensitivity  to  detonation  problem 
requires  a  consideration  of  microscopic  and  macroscopic  characteristics 
of  matter.  The  microscopic  properties  are  described  in  terms  of 
molecular  electronic  structure,  motions  of  nuclei,  and  the  kinds  of 
reactions  which  can  occur.  The  macroscopic  properties  are  related 
to  the  homogeneity  and  heterogeneity  of  the  explosive,  the  mechanical 
response,  which  includes  shear  stresses  and  strains,  and  explosive 
dimensionality  parameters  such  as  critical  diameters  which  give  the 
minimum  charge  diameter  in  which  a  steady  state  detonation  will 
propagate. 

Considerable  progress  has  been  made  in  determining  relationships 
between  sensitivity  and  the  macroscopic  aspects.  It  is  generally 
recognized  that  homogeneous  explosives  initiate  through  a  thermal 
explosion  which  occurs  at  the  entering  face  of  the  acceptor,  where 
the  explosive  has  been  heated  in  its  compressed  state  for  the  longest 
period  of  time.  Substantial  confirmation  of  this  process  has  been 
achieved  by  Hardesty-*-^  in  a  detailed  experimental  study  of  liquid 
nitromethane.  Increasing  the  heterogeneity  of  an  explosive  can 
change  its  sensitivity. 12  For  example,  the  introduction  of  argon 


lOsardo,  R.  D.,  and  Wolfsberg,  M. ,  "The  Wave  Equation  of  a 
Nonlinear  Triatomic  Molecule  and  the  Adiabatic  Correction  to  the 
Born-Oppenheimer  Approximation,"  Journal  of  Chemical  Physics,  Vol.  67 

No.  2.  1977,  p.  593.  - 

l^Hardesty,  D.  R. ,  "An  Investigation  of  the  Shock  Initiation  of 
Liquid  Nitromethane,"  Combustion  and  Flame,  Vol.  27,  1976,  p.  229. 

S.  Army  Materiel  Command,  Engineering  Design  Handbook, 
Principles  of  Explosive  Behavior,  1972,  Chapter  11,  pp,  1-32. 
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gas  into  nitromethane  creates  hot  spots  at  the  positions  of  the 
bubbles,  where  local  ignitions  occur  before  an  entire  cross  section 
of  nitromethane  is  ignited. *2  In  general,  hot  spots,  which  arise 
as  a  consequence  of  shocking,  are  distinguished  as  hydrodynamic 
and  temperature  and  pressure  related,  and  apparently  are  instrumental 
in  causing  heterogeneous  explosives  to  initiate  directly  behind 
the  shock  front,  in  contrast  to  the  situation  found  in  homogeneous 

explosives . 

Bridgman3- 3  has  demonstrated  that  shear  stresses  and  strains  lead 
to  explosive  decomposition  of  hydrostatically  compressed  materials. 

At  the  present  time,  it  is  not  known  exactly  what  role  strain  plays 
in  the  initiation  of  a  shocked  explosive,  although  Teller  has  dealt 
qualitatively  with  this  problem. 1^  The  pressures  in  shocked  materials 
are  generally  believed  to  be  hydrostatic  in  nature,  since  the  yidld 
strengths  are  greatly  exceeded.  These  pressures  relate  to  average 
distributions  of  molecules  which  are  achieved  on  a  time  scale  which 
may  be  long  compared  to  times  associated  with  other  processes  occuring 
on  the  microscopic  scale. 

The  characteristic  times  required  to  yield  commonly  observed 

microscopic  and  macroscopic  effects  in  matter  range  from  10-^6  to 

10“^  sec.  The  orbital  time  for  electronic  motion  is  the  order  of 

10-1°  to  10-3-^  sec.  Since  nuclei  are  much  more  massive  than  electrons 

(m  , ^  /m  .  >  1800)  ,  the  period  for  optic  mode  vibrations  is 

nucleus  electron 

generally  10-3-4  to  10-3-3  sec.,  and  the  period  for  acoustic  mode 
vibrations  is  10-3-3  to  lO”3^  sec.  The  times  required  to  attain  local 
and  global  thermodynamic  equilibria  in  shocked  liquids  are  10-10  sec. 
and  10-7  sec.,  respectively.  The  time  for  nucleation  of  solid 
material  from  the  liquid  state  is  greater  than  the  1-10  ysec.  duration 
in  shock  compression  experiments.  On  the  basis  of  the  above  times, 
then,  any  molecular  model  which  relates  sensitivity  to  molecular 
structure  must  answer  the  question:  why  does  the  initiation  process 
take  times  in  microseconds  as  compared  with  the  much  shorter  times 
for  elementary  motions  within  the  molecules?  Pas tine,  et  al.4  have 
addressed  this  question  in  terms  of  molecular  relaxation  processes. 
Edward  Toton  of  this  Center  is  using  relationships  between  the 
rapid  electronic  motion  on  the  one  hand  and  the  much  slower  nuclear 
motion  on  the  other  to  examine  the  nature  of  the  process  which 


-^Bridgman,  P.  W. ,  "Effects  of  High  Shearing  Stress  Combined 
with  High  Hydrostatic  Pressure,"  Physical  Review,.  Vol.  48,  1935, 
p.  825. 

l4Teller,  E.,  "On  the  Speed  of  Reactions  at  High  Pressures," 
Journal  of  Chemical  Physics,  Vol.  36,  No.  4,  1962,  p.  901. 

^Johnson,  Q.,  Mitchell,  A.,  Keeler,  R.  N. ,  and  Evans,  L. , 

"X-ray  Diffraction  during  Shock  Wave  Compression,"  Physical  Review 
Letters,  Vol.  25,  1970,  p.  1099.  In  this  paper,  the  results  of  the 
X-ray  analysis  of  crystalline  LiF  demonstrate  that  the  uniaxially 
strained  crystal  produced  by  a  130  kbar  shock  relaxes  to  an  isotropic 
structure  in  a  time  which  is  short  compared  to  20  nanoseconds. 
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channels  energy  from  the  lattice  vibrational  modes  of  an  explosive 
into  the  optic  vibrational  modes  of  a  given  molecule.  In  the 
Sections  entitled,  "The  Quantum  Mechanical  Problem"  and  "Calculations 
of  Ground  and  Excited  States  for  the  (N0+)n  Molecular  System"  of  this 
report,  we  shall  begin  to  address  the  above  question  from  the  point 
of  view  of  how  the  energy  within  the  explosive  molecule  is  utilized 
to  break  bonds. 


THE  ACTIVATION  ENERGY  AS  A  POTENTIAL  ENERGY  CONCEPT 


The  reaction  rate  equations  of  a  reacting  fluid  are  given  by 

dt  ZmJ  ^iaij  ,  j  =  1,2, ...N  ,  (1) 

i*l 


where  r  is  the  total  number  of  distinct  reactions  occurring,  a^ .  a 
positive  or  negative  integer  or  rational  number  if  substance  j 
participates  in  the  reaction  i,  and  a^-;  =  0  if  substance  j  does  not 
participate.  Here,  r-  is  a  function  or  the  pressure  P,  temperature 
T,  and  N  composition  variables  n ^ , 

l\  (P,  T,  n^,n2  / . .  •  )  ,  i  —  1,2, ...r  .  (2) 

I\  has  the  general  form 


ri 


ki 


F. 


(3) 


where  k^  is  the  rate  coefficient  for  the  i-th  reaction  and  is  a 
function  of  P  and  T,  anf  is  a  function  of  the  composition 
variables.  The  mathematical  form  of  F.^  contains  products  of  the 
composition  variables.  On  the  other  hand,  k^  describes  the  nature 
of  the  reaction  taking  place  and  utilizes  information  pertaining  to 
molecular  structure  and  the  probabilities  of  intramolecular  and 
intermolecular  energy  transfer.  Since  molecular  structure  and  the 
probabilities  depend  on  P  and  T,  the  fundamental  parameters  contained 
in  k^  of  Equation  (3)  and  the  mathematical  form  of  k^  will  also  vary 
with  P  and  T. 


When  P  and  T  vary  simultaneously,  the  variation  of  Ink.:  with 
1/T  is' 


dlnki 
d  ( 1/T  ) = 


“act,! 


TV 


act,i  dP 
R  dT 


(4) 
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where  R  is  the  gas  constant  and* 


and 


(5) 


(6) 


The  quantities  Eacj;>i  and  Vact>i  are  the  energy  of  activation  and 

the  volume  of  activation,  respectively,  for  the  i-th  reaction.  They 
are  intimately  related  to  molecular  structure.** 

The  analysis  in  Reference  (7)  demonstrates  that  the  function 


E 


int 

act,i 


^act ,  i 


-PV 


act ,  i 


(7) 


is  an  internal  energy  of  activation.  It  is  this  quantity  which  is 
closely  related  to  the  changes  in  potential  energy  when  the  reactants 
pass  from  their  ground  states  to  the  activated  state  in  the  i-th 
reaction. *** 

Integration  of  Equation  (4)  reveals  that  k^  in  Equation  (3) 

contains  exp  (-Eact f i/RT) ,  which  is  the  probability  that  the  reacting 

system  has  an  energy  of  E^}£,i  or  greater.  One  may  then  speak  of 

the  reacting  system  as  having  surmounted  a  barrier  Ea£t,i  in  height. 

This  activation  barrier  is  imbedded  in  a  potential  energy  hypersur¬ 
face  of  dimensionality  given  by  the  geometrical  degrees  of  freedom 
of  the  reacting  system.  The  hypersurface  has  its  origin  in  the 
quantum  mechanical  description  of  the  reacting  system.  Such  a 
description  constitutes  the  most  fundamental  analysis  for  the  system. 
A  discussion  of  the  utility  of  the  potential  energy  hypersurface  in 


*The  variables  V  and  T  may  also  be  chosen.  The  variables  P  and 
T  have  been  selected  here  since  Vact  is  widely  used  by  chemists  to 
interpret  reaction  mechanisms. 

**Since  Eact  is  defined  by  the  partial  derivative  in  Equation  (5)» 
it  is  not  what  one  would  call  the  Arrhenius  activation  energy,  which 
is  defined  by  the  total  derivative  in  Equation  (4) . 

***It  is  noted  here  that  if  P.=  cT  (c  not  a  function  of  T)  is 
substituted  into  Equation  (4) ,  E^£|  becomes  an  Arrhenius  activation 
energy.  However,  P  =  Pq+Pt*  so  that,  in  general,  P  =  cT  will  not 
be  an  adequate  approximation. 
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chemical  kinetics  is  given  by  Laidler.1^*  a  detailed  discussion  of 
these  hypersurfaces  is  presented  in  the  Section  entitled,  "The 
Quantum  Mechanical  Problem." 


THE  QUANTUM  MECHANICAL  PROBLEM 

The  exact  quantum  mechanical  description  of  an  explosive  *  ' 

requires  wave  functions  which  contain  all  of  the  information  about 
the  states  of  the  explosive.  In  a  practical  sense,  however,  it  is 
impossible  and,  frequently,  unnecessary  to  calculate  such  detailed 
information  about  the  system.  Moreover,  even  if  exact  wave  functions 
or  eigenfunctions  of  the  Schrodinger  wave  equation  could  be  obtained, 
they  would  be  so  complicated  that  many  of  the  intuitive  aspects  of 
the  fundamental  interactions  would  be  masked  by  the  quantum  mechanical 
nature  of  the  dynamics.  In  this  report,  models  and  rigorous  computa¬ 
tional  schemes  are  introduced  which  quantitatively  approximate  the 
states  of  highly  compressed  explosives,  and  yet  maintain  a  reasonably 
simple  picture  of  interacting  explosive  molecules. 

Quantum  mechanical  treatments  of  the  fundamental  interactions 
between  N  electrons,  N  nucleons,  N  atoms,  or  N  molecules  have  all 
begun  with  the  averaged  interactions  of  the  N  particles,  where  the 
motion  of  each  particle  is  such  that  it  moves  in  an  averaged  field 
generated  by  the  remaining  N-l  particles.  A  mathematical  basis  for 
this  approach  is  contained  in  the  Hartree-Fock  theory, which  forms 
the  theoretical  framework  for  the  shell  structure  of  nuclei,  atoms, 
and  molecules.  The  shell  model  has  had  immense  success  as  a 
quantitative  and  qualitative  tool  in  the  analysis  of  atoms  and 
molecules.  This  theory  is  as  far  as  one  can  go  in  utilizing  the 
so-called  orbital  picture  in  atoms  and  molecules. 


-*-®Laidler,  K.  J.,  Chemical  Kinetics  (New  York:  McGraw-Hill  Book 
Co.,  Inc.,  1965),  pp.  72-112, 

l7Slater,  J.  C.,  Quantum  Theory  of  Atomic  Structure,  Volume  II, 
(New  York:  McGraw-Hill  Book  Co.,  Inc.,  1960),  pp.  1-30. 

*Frequently,  plots  of  the  free  energy  vs.  reaction  coordinate 
are  used  to  analyze  reaction  paths.  These  plots  represent  the 
totality  of  energy  changes  occurring  during  reaction  and  already 
contain  the  potential  energy  changes  ascribed  to  the  reaction  path 
imbedded  in  the  hypersurface.  The  assumption  of  equilibrium  between 
reactants  and  the  activated  species  in  the  absolute  reaction  rate 
theory  dispenses  with  the  need  to  know  the  pathway  taken  to  reach 
the  activated  state.  Therefore,  when  this  theory  provides  a  good 
description  of  reactions,  one  need  not  be  concerned  with  the  complete 
potential  energy  hypersurface.  In  the  Sections  entitled,  "The 
Quantum  Mechanical  Problem"  and  "Calculations  of  Ground  and  Excited 
States  for  the  (N0+)n  Molecular  System,"  the  surfaces  are  treated  in 
their  entirety. 
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Methods  which  refine  the  Hartree-Fock  theory  have  been  introduced 
which  not  only  improve  the  quantitative  accuracy,  but  which  also 
maintain  intuitive  pictures  of  molecular  and  nuclear  structures.  In 
the  case  of  molecular  electronic  structure,  the  generalized  valence 
bond18  methods,  for  example,  are  finding  increased  usage  in  polyatomic 
molecules.  These  methods  represent  one  aspect  of  the  more  general 
configuration  interaction  techniques  which  treat  the  correlation  of 
motions--o€  electrons  a«d  •of*,  nucleons .  The  generalized  valence  bond 
methods  maintain  the  chemist's  idea  of  the  nature  of  a  covalent 
bond  containing  two  electrons.  The  considerable  quantitative 
improvement  which  these  methods  achieve  is  ascribed  to  the  treatment 
of  the  detailed  interactions  of  any  two  electrons  moving  in  the 
averaged  field  of  the  remaining  N-2  electrons.  Interestingly  enough, 
even  at  the  highest  known  densities  of  matter,  nuclear  theoretical 
techniques  have  demonstrated  that  major  energy  improvements  are  also 
achieved  with  a  pairwise  potential  of  two  nucleons  moving  in  the 
averaged  N-2  nucleon  field.10  Higher-order  corrections  which  provide 
for  the  correlation  of  the  motions  of  electrons  and  of  nucleons  are 
required  to  obtain  complete  agreement  with  experiment;  but,  these 
corrections  obscure  the  traditional  shell  structure  and  electron-pair 
pictures  of  matter. 

The  characterization  of  the  dynamics  of  interacting  molecules 
is  often  not  as  clearcut  as  the  dynamics  of  electrons  and  nucleons, 
since  strongly  interacting  molecules  change  the  electronic  structure 
of  the  individual  molecules.  A  determination  of  equations  of  state 
requires  interaction  potentials.  Jacobs20  has  determined  that  the 
equation  of  state  of  a  strongly  shocked  solid  explosive  is  closer 
to  that  of  a  liquid  than  a  solid.  Whereas  liquids  near  the  triple 
point  are  frequently  analyzed  satisfactorily  in  terms  of  pair 
potentials  4>  ( i ,  j )  between  atoms  and  molecules,  triplet  potentials 
<t>(i,j,k)  and,  perhaps,  higher-order  potentials  are  expected  to  be 
important  at  the  high  densities  of  shocked  explosives  where  the  total 
potential  energy  Uj  for  a  given  configuration  of  atoms  or  molecules 
would  have  the  form 


U 


I 


y]  <f>  (ir  j)  +  2  <Mi»jfk) 
i< j  i< j<k 


+  higher-order 
terms 


(8) 


ISpunning,  Jr.,  T.  H.,  Cartwright,  D,  C,,  Hunt,  W,  J.,  Hay,  P,  J,, 
and  Bobrowicz,  F.  W. ,  "Generalized  Valence  Bond  Calculations  on  the 
Ground  State  (X1Z+Q)  of  Nitrogen,"  Journal  of  Chemical  Physics, 

Vol.  64,  No.  11,  1976,  p.  4755,  and  other  references  cited  therein. 

l0Bethe,  H.  A.,  "Nuclear  Many-Body  Problem,"  Physical  Review, 

Vol.  103,  No.  5,  1956,  p.  1353. 

20Jacobs,  S.  J.,  "On  the  Equation  of  State  for  Detonation 
Products  at  High  Density,"  in  Twelfth  Symposium  (International)  on 
Combustion,  1969,  pp.  501-511, 
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The  summations  in  Equation  (8)  refer  to  the  atoms  or  molecules.  The 
pair  potential  4>'(i,j),  for  example,  would  contain  variables  for  the 
intermolecular  distances  and  relative  orientations  of  two  molecules 
interacting  in  isolation  from  the  other  molecules.  The  actual  pair 
potential  <f>(i,j)  in  the  averaged  field  of  the  remaining  molecules 
would  be  a  simple  function  of  4>'(i,j).  Thus,  a  description  of 
interacting  molecules  can  be  built  up  in  a  manner  analogous  to  that 
done  for  electrons  and  nucleons.  It  will  be  shown  in  the  Section 
entitled,  "Calculations  of  Ground  and  Excited  States  for  the  (N0+)n 
Molecular  System,"  that  the  largest  shifts  in  electronic  excitation 
energies  occur  for  three  interacting  N0+  molecules  so  that  ( i ,  j )  and 
<t>(i>j,k)  are  adequate  for  a  description  of  molecular  electronic 
structure.  The  results  of  detailed  calculations  on  the  system  of 
two  interacting  N0+  molecules  are  presented  in  the  Section  entitled, 
"Calculations  of  Ground  and  Excited  States  for  the  (N0+)n  Molecular 
System,"  from  which  some  idea  of  the  nature  of  cj)(i,j)  is  obtained. 
Triplet  potentials  4>(i,j,k)  will  be  considered  in  a  future  report. 


It  is  frequently  desirable  to  treat  the  interacting  molecules 
as  an  assemblage  of  electrons  and  nuclei  which  yield  a  potential 
energy  Uen.  Once  the  energy  Uen  is  obtained,  it  may  be  used  to 
evaluate  the  pressure  P  from  the  virial  theorem  equation  of  state. 


PV 


(9) 


where  K  is  the  average  total  kinetic  energy,  M  is  the  number  of 
electrons  and  nuclei  in  the  molecular  composite  with  potential 
energy  Uen,  and  "r^  is  a  vector  locating  the  position  of  particle  i. 
Since  we  are  concerned  here  with  electrostatic  forces  among  the 
electrons  and  nuclei,  the  averaged  quantity  in  the  second  term  of 
Equation  (9)  becomes 


M 

-£ 


(10) 


Note  that  the  second  term  of  Equation  (9)  will  be  nonzero  at  0°K  so 
that  P  may  be  written* 


P  =  Pt(V,T)  +  PQ (V) 


(ID 


2lHir schfelder ,  J.  0.,  Curtiss,  C.  F.,  and  Bird,  R.  B., 
Molecular  Theory  of  Gases  and  Liquids,  (New  York:  John  Wiley  & 
Sons,  Inc.,  1954),  pp.  134-135. 

♦Jacobs20  and  Pastine5  have  used  Equation  (11)  ,  which  also 
follows  from  the  Griineisen  equation  of  state. 


13 


NSWC  TR  79-175 


Equations  (9)  and  (10)  may  be  simplified  by  invoking  the  Born- 
Oppenheimer  (BO)  approximation  which  is  permitted  by  the  separation 
of  electronic  and  nuclear  motions.  In  this  approximation,  Slater22 
showed  that 


K 


e 


-E  *i 

nuclei 


(12) 


and 


K  + 
e 


U  =  E 
en  e 


(13) 


where  K  is  the  quantum  mechanical  expectation  value  of  the  electronic 
kineticeenergy ,  and  E  includes  the  total  electronic  energy  and  the 
nuclear  repulsion  enefgy.  Since  K  =  Ke  +  Kn,  substitution  of 
Equations  (12)  and  (13)  into  Equation  (9)  yields 


-  ?E 

nuclei 


(14) 


It  is  seen  here  that  Ee  is  now  the  potential  for  nuclear  motion  to 
be  used  in  the  determination  of  P.  Ee  is  determined  from  the 
electronic  structure  of  the  entire  assemblage  of  interacting  mole¬ 
cules,  as  is  done  in  the  so-called  "supermolecule"  approach  to 
intermolecular  interactions.22  We  now  turn  to  an  analysis  of  the 
nature  of  Ee  and  the  BO  approximation. 

Bardo  and  Wolf sberg'*"®  have  shown  that  the  nonrelativistic 
Schrodinger  wave  equation  in  a  coordinate  system  rotating  with  a 
molecular  system  can  be  written  in  the  form 


(T+V-E)^  + 


j 

s+i 


+  T 


J 

s-1 


*  T+.*s«  +  T--*s-2  ’  0  ■ 


s  =  -J,  -J+l, 


(15) 


22Slater,  J.  C.,  Quantum  Theory  of  Molecules  and  Solids,  Volume  I, 
(New  York:  McGraw-Hill  Book  Co.,  Inc,  1963),  pp.  29-40. 

2 2 Amos ,  A.  T.,  and  Crispin,  R.  J.,  "Calculations  of  Intermolecu¬ 
lar  Interaction  Energies,"  in  Eyring,  H.,  and  Henderson,  D. ,  eds.. 
Theoretical  Chemistry,  Advances  and  Perspectives,  Vol.  2,  (New  York: 
Academic  Press,  1976),  pp.  33-35. 
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In  this  set  of  (2J  +  1)  coupled  equations,  J  and  s  are  the  quantum 
numbers  of  the  total  angular  momentum  and  the  projection  of  the  total 
angular  momentum  onto  the  2-axis  of  a  rotating  coordinate  system, 

respectively.  T,  T+,  T_,  T++,  and  T _  are  kinetic  energy  operators, 

the  forms  of  which  depend  on  the  choice  of  coordinate  system  rotating 

with  the  molecular  system.  T+,  T_,  T++,  T _  arise  from  the  Coriolis- 

like  interactions  between  the  rotational,  electronic,  and  vibrational 
motions.  V  is  the  usual  expression  for  the  potential  energy  which 
includes  the  electrostatic  potential  energy  terms  among  electrons 
and  the  N  nuclei.  E  is  the  total  energy  (excluding  the  energy  of  the 
center-of-mass  motion)  of  the  rotational  and  internal  motions  of  the 
system. 

Equation  (15)  strictly  applies  to  the  entire  assemblage  of 
molecules  in  an  explosive  for  which  J  is  a  good  quantum  number  and, 
in  principle,  provides  an  exact  nonrelativistic  quantum  mechanical 
solution  for  the  explosive  in  any  given  state.  Each  wave  function 
4*3  may  be  expanded  in  terms  of  a  complete  set  of  orthonormal 
electronic  wavefunctions , 


“  X>KS  '  (16 

K 

where  F^g  is  a  function  of  the  3N-6  coordinates  of  the  nuclei,  and 

$K  is  a  function  of  these  nuclear  coordinates  and  the  coordinates  of 
the  electrons.  The  functions  $K  are  orthonormal  in  the  sense  that 

/$*  $_  dv  ..  ,  =  6^.  .  (17 

K  L  electrons  KL 

Substituting  Equation  (16)  into  Equation  (15) ,  multiplying  on  the 
left  by  4>k '  ,  and  integrating  over  ^he  electronic  coordinates,  we  have 


(18) 


<3>k  .  I  T+V  |  $K  ,  >Fr  ,  g  -  EFk,s  +  <‘5K.  IT+I$K.>FK.  >s+1 
+  <$K,  I  T_  I  •  >  Fk  s_i  +  T++Fk  '  ,  s+2  +  T— FK',s-2  = 

“E  I<VIt+vIVfks  +  <V!t+!*k>fk%+i  +  <VlT-1VFK%-ij  . 

K(^K')V 


The  bracket  notation  <| | >  of  Dirac  used  here  denotes  integration  over 
the  electronic  space  and  spin  coordinates  only. 


15 


NSWC  TR  79-175 


The  solution  of  Equation  (18)  is  greatly  simplified  when  the 
functions  are  chosen  in  such  a  way  that  the  number  of  terms  in  the 
summations  in  Equations  (16)  and  (18)  is  as  small  as  possible.  For 
many  physical  problems,  extremely  rapid  convergence  has  been  obtained 
when  the  functions  <t>K  were  chosen  as  electronic  states  within  the 
framework  of  the  BO  approximation,  .  In  this  approximation,  each 
state  ^  of  the  molecular  system  is  very  well  approximated  by 

■  •  • 

*  Tls  * l °  '  <191 


where  L  denotes  additional,  good  quantum  numbers  which  arise  from  the 
approximation  and  which  distinguish  electronic  states,  but  which  are 
nonexistent  in  the  exact  solution  represented  by 

The  equations  of  the  BO  approximation  are  obtained  from  the 
exact  equations,  Equation  (18).  It  was  shown  in  Reference  10  that 
the  operators  T  and  T+  contain  derivatives  with  respect  to  only 
nuclear  coordinates.  In  hartree  units  (1  hartree  =  27.2  ev) ,  the 
operator  T  may  be  written  as 


T  =  T’-I ,  (20) 

n 

where  the  sum  over  n  is  for  the  electrons,  and  is  the  electronic 
Laplacian.  The  operators  T'  ,  T+  and  T+  +  contain  terms  involving 
first  and  second  partial  derivatives,  and  involving  the  reciprocal 
of  the  i-th  nuclear  mass,  One  such  term,  for  example,  is 

(mnAcl,  i)  (32/3x3R|<)  where  R|<  is'an  internuclear  distance  and  x  is  an 
electronic  Cartesian  coordinate.  The  BO  approximation,  then,  utilizes 
the  large  ratio  of  the  nuclear  mass  to  the  electronic  mass,  and  the 
slow  variation  of  the  electronic  wave  function  with  a  nuclear 
displacement  of  less  than  a  bohr  radius  ag  (a0  =  0.529  A.) , 


9<J> 

3R 


BO 

K 


k 


,  BO  | 


(21) 


where  R^  is  an  internuclear  distance  and  9/SR^  is  a  derivative 

M  o  M  0 

contained  in  T ' ,  and  in  T+  =  T+  -  T+ .  T+  and  T+  refer  to  only  the 
nuclear  coordinates  and  the  electronic  coordinates,  respectively. 

T'  and  T^  also  contain  terms  with  electronic  orbital  angular  momentum 

"  -1  DA 

operators  and  mnucl  ^  Consequently,  the  operations  T'4>“  and 
T  +  4>®°  are  small  so  that  Equation  (18)  becomes 
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(T’+EK°-E)FKS  +  ^  Fk!s+1 


+  T 


N 


F 


J 

K ,  S-l 


+  T++  F  K, s+2  +  T~  FK,s-2  =  0 


(22) 


where  we  have  let 


<*k°I-»E^  +  v|C>=ek° 

n 

or 


n 

Equations  (22)  and  (24)  are  the  eigenvalue  equations  of  the  BO 
approximation . 

These  two  equations  represent  a  great  simplification  of  the 
solution  to  the  molecular  problem,  since  the  electronic  motion 
described  by  Equation  (24)  is  treated  independently  of  the  nuclear 
motion  described  by  Equation  (22) .  The  eigenvalues  E^°  are  functions 
of  the  3N-6  nuclear  coordinates  where  N  is  the  number  of  nuclei. 
Equation  (24)  is  solved  for  each  configuration  of  fixed  nuclei. 

Thus,  for  each  such  rigid  configuration,  the  electrons  move  through¬ 
out  the  framework  of  nuclei  and  establish  a  force  field  described 
by  and  E^°  which  acts  on  the  nuclei. Since  the  positions  of 
the  nuclei  are  fixed  in  the  solution  of  Equation  (24),  e|°  is  a 
potential  energy  function  for  nuclear  motion.  This,  of  course,  is 
already  implied  in  the  structure  of  the  operator  T'  +  E®°  in 
Equation  (22),  the  solution  of  which  gives  the  vibrational  motions 
of  the  nuclei  described  by  the  wave  functions  Fj£s,  FK^s±j.f  and 

fk‘?s±2*  For  each  total  energy  eigenvalue  E,  E  -  EE0  is  the 
vibrational  energy. 

Each  function  of  the  nuclear  coordinates  E^°  is  what  is  often 
referred  to  as  an  electronic  state.  The  electronic  state  can 
pertain  to  a  composite  of  molecules  in  the  condensed  phase,  or  to 
a  single  molecule  in  the  gas  phase.  Graphs  of  these  functions 
in  nuclear  coordinate  space  reveal  that  there  is  generally  one 
function  which  gives  energy  values  all  lying  below  the  energy  values 


^4Feynman,  R.  P.,  "Forces  in  Molecules,"  Physical  Review,  Vol  56, 
1939,  p.340. 
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corresponding  to  the  other  functions.  This  lowest  function  is 
ClrileRnthe  ground  electronic  state,  All  higher-lying  functions 

Ey  ,E§  ,  ...  are  called  excited  electronic  states.  The  K-th  graph 
of  the  function  E^  is  a  hypersurface  which  is  a  potential  energy 
surface.  This  is  the  surface  mentioned  in  the  Section  entitled 
"The  Activation  Energy  as  a  Potential  Energy  Concept."  In  the 
case  of  two  interacting  atoms,  as  in  a  diatomic  molecule  or  solid, 
e|°  is  a  curve  (1-dimensional  surface)  which  is  frequently  approxi¬ 
mated  by  a  Morse  potential,  a  Lennard-Jones  (6-12)  potential,  or 
a  Buckingham  potential. 21 

Equation  (24)  may  be  solved  for  a  system  of  reacting  molecules. 
Here  the  function  E®°  gives  a  surface  which  describes  the  reaction 
as  it  proceeds  from  reactants  through  the  activated  state  to  the 
products.  The  reacting  system  is  envisioned  as  a  mathematical 
point  moving  along  some  path  which  originates  at  a  low  point,  the 
position  of  the  reactants,  on  the  surface  and  continues  through 
valleys  between  mountains  and  over  a  pass  into  another  low  point, 
the  position  of  the  products,  on  the  surface.  The  surface  can 
reveal  the  most  probable  course  of  the  reaction  without  having  to 
introduce  the  motions  of  the  nuclei  through  the  solution  of  Equation 
(22).  ,  Unfortunately,  the  generation  of  a  complete  hypersurface  for 
an  electronic  state,  as  it  applies  to  the  total  assemblage  of 
molecules  or  to  even  a  few  reacting  molecules  within  the  assemblage, 
constitutes  an  intractable  computational  problem.  Considerable 
expense  has  been  incurred  in  the  case  of  quantum  mechanical 
calculations  on  three  interacting  hydrogen  molecules.  b  Frequently, 
the  hypersurfaces  for  interacting  and  reacting  molecules  are 
approximated  by  elaborate  potentials  which  contain  Morse  potential 
terms  and  exponential  repulsion  terms  with  coefficients  which  modify 
the  attractive  parts  of  Morse  potentials  so  that  exothermicity  in 
reactions  is  achieved.  These  potentials  have  been  used  in  conjunction 
with  Monte-Carlo  calculations  on  reaction  trajectories.  6 

It  is  often  unnecessary  to  consider  hypersurfaces  in  their 
entirety.  For  example,  the  analysis  of  the  breaking  of  a  single 
bond  in  a  molecule  may  be  carried  out  by  examining  a  cut-out  of  the 
surface  which  corresponds  to  the  variation  of  a  single  bond  distance. 
In  the  case  of  the  sensitivity  of  an  explosive,  the  bond-breaking 
process  may  be  in  some  cases  the  most  important  step  in  the  reaction 
mechanism  for  initiation  of  detonation.  It  is  this  process  which  is 
to  be  examined  in  detail  in  this  report. 


25Ross,  M. ,  and  Shishkevish,  C.,  "Molecular  and  Metallic 
Hydrogen,"  R-2056-ARPA,  1977. 

26Blais,  N.  C.,  and  Bunker,  D.  L.,  "Monte  Carlo  Calculations.il. 
The  Reactions  of  Alkali  Atoms  with  Methyl  Iodide,"  Journal  of  Chemical 
Physics,  Vol.  37,  No.  11,  1962,  p.  2713.  See  also  Karo,  A.  M. ,  and 
Hardy,  J.  R. ,  "Molecular  Dynamics  of  Shock-Initiated  Detonations," 
International  Journal  of  Quantum  Chemistry,  Vol.  12,  Suppl.  1,  1977, 
p.  333. 
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The  BO  approximation  is  really  a  semiclassical  treatment  made 
possible  by  the  large  nuclear- to-electron  mass  ratio.  The  solution 
of  Equation  (24)  constitutes  a  violation  of  the  quantum  mechanical 
nature  of  motions  of  nuclei,  since  the  nuclei  are  held  in  fixed 
positions,  giving  each  nucleus  a  simultaneously  known  position  and 
velocity  (velocity  =  0)  which  violates  the  uncertainty  principle. 

It  is  this  aspect  which  permits  the  use  of  Monte-Carlo  calculations 
in  the  analysis  of  reaction  mechanisms.  As  long  as  single  surfaces 
E^°  are  sufficient  for  describing  reactions,  the  BO  approximation 
is  valid  to  a  high  degree  of  approximation.  The  semiclassical 
treatment  is  less  accurate  when  the  molecules  contain  nuclei  of 
small  mass.  Such  is  the  case  for  isotopic  exchange  reactions 
involving  the  hydrides  and  deuterides,  where  quantum  mechanical 
effects  due  to  the  interactions  of  nuclear  and  electronic  motions 
become  important.  Here,  the  interactions  are  contained  in  a 
first-order  correction  term 


'K 


=  <$ 


BO  | 
K  ! 


$B0> 

K  ' 


(25) 


which  was  dropped  in  Equation  (22)  ,  but  which  adds  on  to  E®°  when  the 
BO  approximation  no  longer  holds.  Since  the  BO  approximation  is 
frequently  referred  to  as  an  adiabatic  theory,  CK  is  called  the 
adiabatic  correction,  since  a  new  potential  energy  surface  E^°  +  Cr 
is  considered.  ^7 

So  far  we  have  considered  only  those  physical  situations  where 
the  surfaces  e|°  are  sufficiently  far  apart  so  that  only  one  surface 
is  needed  in  the  analysis  of  a  given  reaction.  It  is  frequently 
found,  however,  that  the  shapes  of  the  surfaces  vary  enough  with 
nuclear  configuration  so  that  two  or  more  surfaces  come  close  to 
intersecting  or,  in  fact,  do  intersect.  The  situation  where  two 
curves  closely  approach  each  other  without  actually  crossing  is 
often  referred  to  as  an  "avoided  crossing."  If  two  surfaces  intersect 
at  a  point  corresponding  to  a  particular  configuration  of  the  nuclei, 
we  have  the  condition 


(26) 


27sardo,  R.  D. ,  and  Wolfsberg,  M. ,  "A  Theoretical  Calculation  of 
the  Equilibrium  Constant  for  the  Isotopic  Exchange  Reaction  between 
H2O  and  HD,"  Journal  of  Physical  Chemistry,  Vol.  80,  No.  10,  1976, 
p.  1068.  It  was  shown  in  this  paper  that  the  adiabatic  correction 
leads  to  equilibrium  constant®  which  differ  from  the  experimental 
ones  by  less  than  the  experimental  error. 
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Teller2®  has  shown  that  the  shapes  of  two  surfaces  in  the  neighborhood 
of  an  intersection  define  a  double  cone  when  each  surface  is  a 
function  of  two  parameters.  More  complicated  types  of  intersections 
are  discussed  by  Carrington. These  intersections  can  have  a 
profound  effect  on  reactions  since  the  BO  approximation  breaks  down 
and  Equations  (21),  (22),  and  (24)  no  longer  apply.  Here  we  find 
that  the  semiclassical  approach  fails,  and  that  the  quantum  mechanical 
effects  manifest  themselves  completely  in  Equation  (18) . 

The  theory  of  radiationless  transitions2®  is  based  on  the 
intersections  of  surfaces.  Intersections  among  the  excited  state 
surfaces  are  much  more  prevalent  than  among  these  surfaces  and  the 
ground  state  surface.  Hence,  radiationless  processes  (e.g., 
predissociation)  and  photochemical  reactions  frequently  involve 
excited  electronic  states.  For  these  reactions,  electronic  excitation 
is  first  required  to  reach  an  excited  state  from  which  the  radiation¬ 
less  transition  occurs.  It  will  be  demonstrated  later  that  surfaces 
belonging  to  excited  electronic  states  intersect  ground  state  surfaces 
during  compression  and,  therefore,  may  affect  the  nature  of  the 
bond-breaking  reaction  during  initiation. 

The  separation  between  two  surfaces  E^°  and  E§°  is  given  by 


AE  =  <Y(H11-H22)2  +  4H22  ,  (27) 

where  H^j  is  the  matrix  element 

Hij  =  <Vt'iiHe|4,j>  '  (28) 


He  is  the  electronic  hamiltonian  of  Equation  (24) .  The  eigenfunctions 

and  iJ>^0  belong  to  the  symmetry  group  of  He.  The  wave  functions 

¥i  and  are  basis  functions.  Longuet-Higgins2-*-  has  proven  that  the 
necessary  and  sufficient  conditions  for  an  intersection,  AE  =  0 ,  of 
two  surfaces  to  occur  are,  from  Equation  (27), 


28 

Teller,  E.,  "The  Crossing  of  Potential  Surfaces,"  Journal  of 
Physical  Chemistry,  Vol.  41,  No.  1,  1937,  p.  109. 

2yCarrington,  T.,  "The  Geometry  of  Intersecting  Potential 
Surfaces,"  Accounts  of  Chemical  Research,  Vol.  7,  1974,  p.  20. 

2®Jortner,  J.,  Rice,  S.  A.,  and  Hochstrasser ,  R.  M. , 
"Radiationless  Transitions  in  Photochemistry,"  in  Noyes,  W.  A.  et  al. , 
eds.,  Advances  in  Photochemistry,  Vol.  7,  (New  York:  Wiley,  1969), 
pp.  149-309. 

21Longuet-Higgins ,  H.  C.,  "The  Intersection  of  Potential  Energy 
Surfaces  in  Polyatomic  Molecules,"  Proceedings  of  the  Royal  Society, 
London ,  Vol.  A344,  1975,  p.  147. 
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H11  =  H2 2  '  Hi2  =  0  •  (29) 

In  general, 


H11_H22  =  G  '  H12  =  G’  '  (30) 

where  both  G  and  G’  are  functions  of  all  of  the  nuclear  coordinates. 

If  these  coordinates  are  internuclear  distances  and  one  distance  is 
varied  while  all  others  are  fixed,  a  curve  is  generated  which  is 
imbedded  in  the  surface  under  consideration.  Von  Neumann  and  Wigner32 
introduced  the  noncrossing  rule  for  two  such  curves  when  they  showed 
that  H12  5*  0  when  the  basis  functions  H'i  and  have  the  same 
symmetry  (spin,  parity,  angular  momentum)  so  that  AE  never  vanishes. 

In  the  case  of  a  diatomic  molecule,  there  is  only  one  internuclear 
distance  parameter  for  which  it  is  impossible  to  satisfy  both 
conditions  in  Equation  (29)  by  solving  the  two  equations,  Equation 
(30),  simultaneously.33  jn  the  case  of  a  polyatomic  molecule  with 
two  or  more  independent  internuclear  distance  parameters,  at  least 
two  of  these  parameters  may  be  found  which  satisfy  the  two  conditions 
of  Equation  (29)  simultaneously.  If  H^2  is  complex  when,  for  example, 
magnetic  forces  are  important,  additional  parameters  are  introduced. 
Hence,  at  least  as  far  as  the  geometrical  parameters  are  concerned, 
intersections  between  two  surfaces  corresponding  to  electronic  states 
of  the  same  symmetry  may  occur  for  polyatomic  molecules,  as  opposed 
to  the  situation  for  diatomics.  The  major  point  here  is  that  there 
are  no  symmetry  restrictions  in  polyatomic  molecules  which  would 
prevent  the  intersections  of  surfaces,  so  that  one  may  expect 
intersections  to  be  prevalent  in  large  molecules.  When  diatomic 
molecules  interact  with  each  other  at  high  compressions,  one  must 
really  consider  the  system  to  be  a  "supermolecule"  which  is, 
effectively,  a  "polyatomic  molecule,"  The  noncrossing  rule  may 
break  down  for  diatomics  at  high  compressions. 


32von  Neumann,  J.  and  Wigner,  E.,  "Uber  das  Verhalten  von 
Eigenwerten  bei  adiabatischen  Prozessen,"  Zeitschrift  fur  Physik, 

Vol.  30,  1929,  p.  467.  ' 

33Hatton,  G.  J.,  "The  Noncrossing  Rule  and  Spurious  Avoided 
Crossings,"  Physical  Review  A,  Vol.  14,  No.  3,  1976,  p.  901.  Hatton 
has  suggested  that  there  may  be  dynamical  symmetries  besides  the 
spatial  and  spin  symmetries  which  establish  additional  relationships 
among  matrix  elements,  even  though  the  states  involved  have  the  same 
spatial  and  spin  symmetries.  Here,  the  traditional  noncrossing  rule 
of  von  Neumann  and  Wigner  would  no  longer  be  applicable.  The  condition 
H12  =  0  exists  for  all  symmetries  as  a  consequence  of  the  following 
important  theorem:  If  there  exists  a  complete  set  of  commuting 
operators  which  includes  He,  and  4*$.  and  are  eigenfunctions  of  at 
least  one  operator  (not  He) ,  then  H12  =  0  unless  the  two  eigenvalues 
a^  and  aj  satisfy  a^  =  aj. 
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Equation  (21)  is  no  longer  satisfied  when  the  eigenfunctions 
of  H~  in  Equation  (24)  vary  rapidly  near  an  "avoided  crossing." 
Here,  tne  designation  BO  has  been  dropped  since  it  no  longer  applies 
to  <fK  and  E^  when  the  matrix  elements  * |t 1 +T++T_ | $K >  in  Equation 


(18)  are  no  longer  negligible.  However,  a  new  set  of  basis  functions 
l*K\  may  be  found  which  are  different  from  and  which  satisfy 


(31) 


The  wave  functions  are  no  longer  eigenfunctions  of  He,  but  still 
allow  a  separation  of  motions  of  nuclei  and  electrons.  In  this  case, 
Equation  (18)  becomes 


(  T'+4,k,-E)F^,s  +  T^,;S+1  + 


N=J 


-J 


+  T  F^  +  T 

i++t K 1 , s+2  --  K' ,s-2 


--E 

K(^K’) 


FJ 

K  Ks 


(32) 


where 


^K'  ^K'K 


(33) 


It  is  noted  here  that 


♦k  -  Ev,  -  (34) 

L 

where  the  coefficients  UKL  comprise  a  unitary  transformation,  and 
are  functions  of  the  nuclear  coordinates.  When  the  nuclear 
configurations  are  far  removed  from  the  one  at  the  "avoided  crossing," 

UKL  ■  5KL- 

Diagrams  of  an  "avoided  crossing"  for  two  curves  are  presented 
in  Figure  la.  Figure  la  shows  thfe  situation  when  the  tfco  curves 
Ei  and  E2  closely  approach  each  other  and  bend  sharply  away.  The 
derivative  in  Equation  (21)  is  large  in_the  region  of  closest  approach. 
In  Figure  lb,  the  new  choice  of  and  $2  smooths  out  the  curves  and 
leads  to  an  intersection,  =  £22*  Equation  (31)  is  obeyed  in 

Figure  lb. 
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It  is  apparent  from  Equations  (32)  and  (33)  that  the  potential 
energy  "seen"  by  the  nuclei  has  meaning  only  when  the  kinetic  energy 
of  the  nuclei  is  considered  and,  therefore,  only  with  respect  to  the 
total  energy  E  attained.  If  E  has  values  close  to  <£kK (=£ll) • 
molecular  system  has  a  good  chance  of  "escaping"  through  an  intersec¬ 
tion  and  of  achieving  nuclear  configurations  beyond  the  intersection 
by  moving  on  E^.  Here,  the  system  may  have  sufficient  time  to 
adjust  even  though  |34>k/8R]c|  may  be  large.  On  the  other  hand,  if  E 
has  values  close  to  E^(=El),  which  is  an  intersection  for  the 
eigenvalues  of  Equation  (24) ,  the  system  has  a  good  chance  of 
"escaping"  on  <fKK  to  nuclear  configurations  beyond  the  intersection 

even  though  | Rk I  maY  large.  When  E  is  much  larger  or  much 

smaller  than  £kk=£ll'  the  surface  "selected"  for  the  most  probable 
motion  is  either  Er  or  El« 

The  concept  of  the  probability  of  "surface  hopping,"  which  was 
dexcribed  in  the  previous  paragraph,  was  introduced  by  Zener,  and 
Landau, and  Stiickelberg , 3 6  and  is  discussed  by  Eyring,  et  al.33 
within  the  context  of  the  general  theory  of  reaction  rates.  These 
authors  consider  two  surfaces  given  by  E^  and  E2  of  the  molecular 
system  "hopping"  from  E^  to  E2/  E^<E2.  Using  the  wave  functions 
and  energies  defined  previously  in  this  report,  we  have  at  =  £22  > 


♦1.2-  ^(*1**2) 

and 


(35) 


He$l  ^11*1  "  ^12^2 


(36) 


He^2  ^22^2  "  ^12^1 


(37) 


where  is  given  by  Equation  (33).  The  expression  for  the 

probability  of  "surface  hopping,"  the  transition  probability,  is 


34zener,  C. ,  "Non-Adiabatic  Crossing  of  Energy  Levels," 
Proceedings  of  the  Royal  Society  of  London,  Ser.  A,  Vol.  137,  1932, 
P.  696.  ~  ‘  "  . . . 

-^Landau,  L.  D.,  and  Lifshitz,  E.  M. ,  Quantum  Mechanics ,  2nd 
Edition  (Reading,  Mass.?  Addison-Wesley  Publishing  Co. ,  Inc. ,  1965), 
no.  322-330. 

36Stuckelberg,  E.  C.  G.,  "Theorie  der  Unelastischen  Stosse 
Zwischen  Atomen,"  Helvetica  Physica  Acta,  Vol.  5,  1932,  p.  369. 

33Eyring,  H.,  Walter,  J. ,  and  Kimball,  G.  E.,  Quantum  Chemistry 
(New  York:  John  Wiley  &  Sons,  Inc.,  1944),  pp.  299-331 . 
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P  =  exp 


«"2*12 

hv|S;L-s2|  J 


(38) 


where  h  is  Planck's  constant,  v  is  the  nuclear  velocity  along  a  curve, 
and  s^  =  S2  =  3<^22/^Rk  are  fcbe  slopes  evaluated  at  the 

nuclear  configuration  where  =  <?22*  T*le  probability  that  the 

system  stays  on  is 


P '  =  1-P 


(39) 


Equations  (38)  and  (39)  tell  us  that  for  a  high  nuclear  velocity  and 
a  small  separation  of  curves,  1 2  £^2 1  ,  at.  the  "avoided  crossing,"  the 
probability  is  large  that  the  system  "hops,"  and  that  P'  is  small. 
Equations  (38)  and  (39)  continue  to  apply  when  the  surfaces  E^  and 
E2  themselves  intersect  when,  for  example,  $1  and  $2  are  of  different 
symmetry. 

The  preceding  discussion  of  "avoided  crossings"  was  concerned 
with  surfaces  Ek  and  EK 1  which  are  the  eigenvalues  in  Equation  (24). 
The  corresponding  eigenfunctions  $K  and  $K 1  have  the  same  symmetry 
so  that  the  noncrossing  rule  is  obeyed.  Equations  (31)  -  (39)  arose 
as  a  consequence  of  the  system  being  able  to  "hop"  from  Er  to  ER , . 

In  the  case  of  an  intersection  Ej<  =  EK 1  ,  <J>K  and  Or  <  will  usually 
satisfy  Equation  (21) ,  but  not  Equation  (22) .  At  the  intersection 
Eg  =  E2 /  degenerate  perturbation  theory  is  used  which  involves 
diagonalizing  a  two-dimensional  matrix  with  matrix  elements 

|  T  ' +T++T_  |  4>  j  y.  This  procedure  yields  two  new  wave  functions 


ii$i  +  ui 2^  2 

t 

(40) 

21$1  +  U22$2 

• 

(41) 

These  wave  functions  are  associated  with  two  new  surfaces  which  do 
not  intersect.  It  should  be  pointed  out  that  the  unitary  transforma¬ 
tion  of  Equations  (40)  and  (41)  is,  in  general,  different  from  that 
defined  in  Equation  (34).  The  transformation  in  Equation  (34) 
defines  the  wave  functions  $1  and  $2  which  minimize  <$1 1 T '+T++T_ | $2 > 
so  that  Equation  (32)  is  satisfied.  Diagonalization  of  the  matrix 
itself  does  not  guarantee  the  smallness  of  matrix  elements  in  the 
different  basis  set.  A  complete  analysis  of  a  given  physical  problem 
requires  the  use  of  the  exact  wave  equations.  Equation  (18),  in  a 
perturbation  treatment.  Once  a  given  set  of  surfaces  is  specified, 
Equations  (38)  and  (39)  are  applicable.  In  general,  the  matrix 
element  Equation  (38)  is  given  by 
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<1^  |t'+t++tJ  y2> 


(42) 


and  not  Equation  (33).  In  a  future  report,  we  consider  the 
mathematical  development  and  computation  of  £l2* 

The  importance  of  the  preceding  discussion  for  the  shock 
initiation  process  is  that  compression  shifts  and  warps  surfaces 
and,  consequently,  may  create  intersections  among  the  excited  state 
surfaces  E]_,  E2,  and  the  ground  state  surface  Eo?  or  bring  the 

surfaces  into  close  proximity.  New  channels  for  reaction  may  be 
created  in  the  process.  As  a  consequence,  alterations  of  the 
surfaces  will  be  manifested  in  changes  in  reaction  kinetics  and  in 
frequency  and  bandwidth  shifts  for  infrared,  visible,  and  ultraviolet 
spectra.  If  by  excited  and  ground  state  surfaces  we  mean  those  of 
the  isolated  molecule,  then  the  compression  shifts  and  warps  these 
surfaces  in  a  perturbative  sense  as  the  molecules  approach  each 
other.  If  we  mean  that  the  surfaces  belong  to  a  composite  system 
of  molecules,  the  surfaces  for  the  composite  will  also  be  shifted 
and  warped  to  an  extent  depending  on  the  state  of  compression.  At 
high  compressions,  one  can  really  consider  only  the  surfaces  for  the 
composite  system.  In  this  case,  a  perturbation  analysis  yields 
either  slow  convergence  or  divergence.  However,  an  "effective  or 
average  set  of  surfaces  per  molecule"  may  be  defined  by  dividing 
the  energies  Ej(  by  the  number  of  molecules  in  the  composite. 

The  shifting  of  existing  intersections  and  the  creation  of  new 
intersections  will  change  activation  energies  from  their  values 
found  under  ambient  conditions.  These  activation  energies  correspond 
to  the  ones  discussed  in  the  Section  entitled  "The  Activation  Energy 
as  a  Potential  Energy  Concept"  and  designated  there  as  E|2£*j_  in 
Equation  (7).  The  solutions  of  Equation  (24),  which  yield  surfaces 
at  various  stages  of  compression  of  a  molecular  system,  can  be  used 
to  compute  activation  energies  which  are  accurate  to  within  10%  of 
the  real  values.  For  example,  a  possible  activation  energy  is 

Ea"t,i  “E0(x>  ’  V0)  -  <43> 

where  Eq(X)  and  Eg(0)  are  the  eigenvalues  of  Equation  (24)  at  the 
nuclear  configurations  corresponding  to  the  activated  system  at  the 
intersection  and  to  the  reactants  at  equilibrium,  respectively.  In 
order  to  obtain  elementary  activation  energies  to  a  higher  degree  of 
accuracy,  and  to  complete  the  description  of  the  reaction  kinetics, 
higher-order  contributions  must  be  generated  from  Equation  (18) . 

These  contributions  would  yield  c^2  (Equation  (42))  for  use  in 
Equations  (38)  and  (39).  Then,  when  the  rate  coefficient  ki  is 
expressed  as 
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(44) 


we  have  a  description  of  the  reaction  in  terms  of  the  activation 
energy  E^{3|  j ,  and  of  the  probability  P',  Equation  (39),  averaged 
over  the  velocities  v.  E3jg|  ^  in  Equation  (44)  is  the  activation 
energy  per  molecule.  In  the  case  of  an  intersection  Eg  *  E]_, 
Equation  (44)  indicates  that,  for  a  given  nuclear  velocity  v  and 
a  slope  difference  | s S2 1 ,  an  increase  in  £12  and  a  corresponding 
decrease  in  Ej{3|  i  greatly  increases  reaction  rates.  This  is  the 
behavior  which  may  yield  the  fast  reaction^  in  detonations. 

It  is  seen  that  the  analysis  consists  of  two  parts.  First, 
Equation  (24)  is  solved  which  yields  Ef(  at  various  compressions. 

The  existence  of  an  intersection  Ej<  =  Ejj  1  permits  the  determination 
of  a  good  approximation  to  the  activation  energy  given  by  Equation 
(43).  Second,  the  utilization  of  perturbation  theory  in  conjunction 
with  Equation  (42)  refines  the  picture  in  terms  of  Equation  (44) . 

In  the  next  section,  the  effects  of  compression  on  EK  are  analyzed, 
and  E^2t,i  *-s  evaluated  for  the  system  of  interacting  N0+  molecules. 

CALCULATIONS  OF  GROUND  AND  EXCITED  STATES  FOR  THE  (N0+) n  MOLECULAR 
SYSTEM 


The  solutions  to  Equation  (24)  are  obtained  by  configuration 
interaction  (Cl)  methods.  Such  solutions  are  of  much  higher  accuracy 
than  is  often  needed  in  calculations,  and  are  obtained  at  great 
computer  expense.  Frequently,  one  looks  for  techniques  which  reveal 
as  much  information  as  is  required  for  lower  cost.  Some  of  these 
techniques  are  briefly  discussed  in  Appendix  A.  Instead  of  solving 
Equation  (24)  ,  these  methods  yield  solutions  to  the  new  eigenvalue 
equation 


(He-V+VaPP^K  =  Vk 


(45) 
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where  Vapp  is  an  approximate  potential  energy  operator.  Vapp  gives 

the  interaction  of  single  electrons  with  averaged  fields  in  the  case 
of  the  LCAO-MO-SCF  method38  or  the  interaction  of  electrons  which 
modify  the  averaged  fields  in  the  case  of  the  CNDO/S-CI  method.39 
The  meanings  of  the  acronyms  are  given  in  Appendix  A.  Calculations 
are  now  discussed  which  utilize  these  methods  to  provide  solutions 
to  Equation  (45)  at  relatively  low  cost. 

The  molecular  system  to  be  considered  in  this  report  is  (N0+)n 
where  N0+  is  the  closed-shell  counterpart  of  the  explosive  open-shell 
NO  molecule,  and  n  is  the  number  of  molecules  required  to  simulate 
compressed  states.  The  system  (N0+)  n  is  really  a  supermolecule  on 
which  LCAO-MO-SCF  and  CNDO/S-CI  calculations  are  carried  out. 

The  (N0+)n  molecular  system  represents  a  compromise  between  the 
NO  explosive  and  the  more  common  explosive  molecular  system.  A 
description  of  the  NO  liquid  explosive  is  complicated  by  dimeriza¬ 
tion.4^  However,  most  explosive  molecules  have  closed  shells  and, 
therefore,  do  not  dimerize  under  ambient  conditions.  The  more 
complicated  theoretical  problem  concerning  (N0)n  will  be  considered 
in  a  future  report. 

The  CNDO/S-CI  method  was  used  on  (N0+)n,  n  =  1,2, 3, 4, 5,  and  6, 
in  order  to  ascertain  the  number  of  molecules  required  to  yield  the 
largest  shifts  in  electronic  excitation  energies.  Various  configur¬ 
ations  of  N0+  molecules  were  considered  in  each  grouping.  The 
results  of  the  calculations  are  displayed  in  Figure  2  where  the 
excitation  energies  L&  from  the  ground  state  to  the  lowest  excited 
state  are  plotted  against  n.  The  graphs  represent  the  variation  of 
A  £  in  kcal  per  mole  of  supermolecules  (N0+)n.  In  all  calculations, 
the  N-0  bond  distance  in  N0+  was  fixed  at  its  equilibrium  gas  phase 
value  of  1.062&.  The  straight  line  pertains  to  all  N0+  molecules 
separated  by  distances  greater  than  3%.  between  nearest  atoms, 
regardless  of  molecular  orientation.  The  3&  distance  is  a  lower 
bound  to  that  found  in  the  condensed  system  at  ambient  pressures. 

This  line  proves  that  under  normal  conditions  the  condensed  phase 
spectra  are  almost  identical  to  gas  phase  spectra.  The  remaining 
curves  indicate  the  variability  of  a£  for  N0+  molecules  in  linear 


38Roothaan,  C.  C.  J. ,  "New  Developments  in  Molecular  Orbital 
Theory."  Reviews  of  Modern  Physics,  Vol.  23,  No.  2,  1951,  p.  69. 

39Del  Bene,  J.,  and  Jaffe,  H.  H.,  "Use  of  the  CNDO  Method  in 
Spectroscopy.  I.  Benzene,  Pyridine,  and  the  Diazines,"  Journal  of 
Chemical  Physics,  Vol.  48,  No.  4,  1968,  p.  1807,  and  references 
cited  therein. 

480hlsen,  J.  R. ,  and  Laane,  J.,  "Characterization  of  the 
Asymmetric  Nitric  Oxide  Dimer  0=N-0=Nby  Resonance  Raman  and 
Infrared  Spectroscopy,"  Journal  of  the  American  Physical  Society, 
Vol.  100,  No.  22,  1978,  p.  6948. 
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and  zigzag  chains  where  the  nearest  atoms  of  neighboring  molecules 
are  separated  by  R  =  2A.  The  nearest  neighbor  distance  for  the 
parallel  system  is  R  =  2.3&.'  In  the  linear  and  zigzag  chains,  the 
R  =  2&.  separation  corresponds  to  distances  of  nearly  3&  between  the 
centers-of-mass  of  the  individual  molecules.  It  will  be  shown  later 
that  this  separation  corresponds  to  a  pressure  P0  a  136  kbar  on  the 
0°K  isotherm.  In  the  case  of  parallel  molecules,  R  =  2.3A  is  also 
the  distance  between  centers-of-mass. 

It  is  interesting  to  note  that  a£“  reaches  a  minimum  for  three 
molecules  at  R  =  2 A,  regardless  of  orientation,  and  approaches  a 
constant  value  as  n  +  «>.  The  reason  for  this  behavior  is  manifested 
in  the  nodal  and  antinodal  (bonding)  structure  of  the  excited  state 
molecular  orbitals.  Each  of  the  molecular  orbitals  for  the  electronic 
states  has  m  nodes  which  oppose  bonding,  and  m'  partial  bonds  between 
molecules.  The  number  of  nodes  m  is  nearly  proportional  to  the 
number  of  N0+  molecules,  and  m*  is  nearly  proportional  to  the  number 
of  pairs  of  N0+  molecules.  The  ratio  r  =  m/m'  is  constant  for  the 
molecular  orbitals  in  the  ground  state  wave  function.  On  the  other 
hand,  r  for  the  excited  state  tt*  molecular  orbital,  for  example, 
reaches  a  minimum  of  3/2  for  (N0+)3,  and  is  larger  for  n  /  3.  In 
all  cases,  the  nodal  integrity  of  each  N0+  "molecule"  is  maintained. 

The  constant  value  of  A<?  for  large  n  lies  close  to  the  value  of 
A<£  =  125  kcal/mole  for  (N0+)2-  This  result  suggests  that  the  shifts 
of  excited  states  for  an  explosive  are  primarily  functions  of  the 
interactions  of  two  or,  perhaps,  three  molecules.  This  feature 
implies  that  extensive  calculations  of  excited  states  may  be 
performed  on  only  a  small  number  of  molecules.  These  results  are 
expressed  even  more  clearly  in  Figure  3  where  h€,  given  in  kcal 
per  mole  of  NO*  molecules,  is  seen  to  approach  constant  values  in 
accordance  with  Figure  2.  Long-range  effects  are  seen  to  be  relatively 
unimportant  for  all  structures.  These  effects  tend  to  be  greatest  for 
the  linear  chain  since  electron  mobility  among  the  molecules  is 
largest  here. 

In  accordance  with  the  BO  approximation  and  the  corresponding 
Franck-Condon  principle,  the  excitation  energies  a£  were  computed 
at  fixed  nuclear  configurations.  Unfortunately,  limitations  within 
the  CNDO/S-CI  method  preclude  the  computation  of  surfaces  when  the 
bond  distance  in  N0+  is  substantially  greater  than  its  equilibrium 
value  of  1.062^.  We  now  describe  surfaces  generated  by  the  LCAO-MO- 
SCF  method. 

As  mentioned  before,  the  initiation  of  detonation  process 
involves  the  breaking  of  bonds.  We  are  interested  in  the  mechanism 
which  accelerates  the  rate  of  bond  scission.  To  this  end,  the 
LCAO-MO-SCF  method  was  used  to  generate  cuts  out  of  the  total 
reaction  surface  of  six  dimensions  for  (NO+)2»  One  of  these  cuts, 
a  curve,  corresponds  to  the  change  in  potential  energy  as  N-0  bonds 
are  stretched  while  the  other  internuclear  distances  are  held  fixed. 

To  a  good  first  approximation,  the  (NO+) 2  system  should  reflect 
changes  in  the  N-0  bond  strength,  since  Figures  2  and  3  show  that 
this  system  yields  the  largest  depression  of  from  the  single  NO+ 
molecular  value  of  173  kcal/mole. 
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In  Figure  4,  the  ground  state  potential  energy  curve  c  q  (1G) 
for  the  isolated  molecule  is  plotted  against  the  N-0  bond  distance 
R'  expressed  in  bohr  radii  )1  bohr  =  0.529A).  All  energies  are 
referred  to  the  minimum  at  R^j^  =  2.007  bohrs  (1.062A).  The  section 
of  the  curve  for  R  <1^10  corresponds  to  the  repulsion  between  atoms. 
Extension  of  the  curve  would  show  that  -*■  00 .  Here,  the  repulsion 

of  the  nuclei  dominates  ^q*0.  The  section  of  the  curve  for  R  <R  min 
shows  a  steep  rise  and  an  eventual  falling  off  to  a  limit  character¬ 
istic  of  the  energy  of  the  atoms  N  and  0+.  The  shape  of  the  curve 
is  dictated  bv  the  virial  theorem  relationship  of  Equation  (12) 
where  Ee  =  @  o°  •  The  distance  from  the  bottom  of  the  potential  well 

at  Rmin  to  the  plateau  f<?r  dissociation  is  approximately  254  kcal/mole. 
Thisvalue  pertains  to  E^ct  an<^  slmHar  to  the  bond  dissociation 
energy,  256  kcal/mole,  or  CO,  which  is  isoelectronic  with  NO+. 


The  two  lowest  excited  state  curves  ^^(lE)  and  <?^°(1E')  for 
NO  are  also  plotted  in  Figure  4.  These  curves  were  obtained  from 
Equation  (A7)  in  Appendix  A,  using  the  MO 1 s  from  the  ground  state 
calculation.  Since  additional  computer  programming  is  required,  no 
SCF  calculations  on  the  excited  states  were  performed  here.  An 
SCF  calculation  could  lower  each  curve  by  as  much  as  30  kcal/mole. 

It  is  noted  from  the  shapes  of  the  curves  that  NO+  is  bound  in  these 
excited  states.  The  dissociation  limits  of  and  were  not 

computed. 


Figure  5  gives  the  orientations  of  the  two  NO+  molecules  in 
(NO+)2  for  a  head-on  collision  and  a  broadside  collision.  These  two 
configurations  do  not  necessarily  give  the  lowest  energies  for  the 
(NO+) 2  system.  They  are  considered  here  for  direct  comparison  with 
the  (NO) 2  system  which  is  a  dimer  with  the  indicated  linear 
configuration  (head-on  collision)  at  ambient  pressures. 40  This 
system  will  be  discussed  in  a  future  report.  The  broadside 
collision  provides  a  direct  reaction  path  for  generation  of  the 
detonation  reaction  products  N2  and  O2  in  the  NO  explosive. 41  An 
LCAO-MO-SCF  calculation  was  performed  at  each  set  of  values  for 
R  and  R1  where  R  and  R'  were  chosen  in  the  ranges  1  <  R  £  4& 

(1.9  <  R  *  7.6  bohrs)  and  0.9  £  R’  £  1.6&  (1.7  <  R*  £  3.0  bohrs). 

The  resulting  curves  are  plotted  in  Figures  6  and  7. 


4lRamsay,  J.  B.  ,  and  Chiles,  W.  C.,  ’’Detonation  Characteristics 
of  Liquid  Nitric  Oxide,"  in  Proceedings  of  the  Sixth  Symposium 
(International)  on  Detonation")  24  August  1976,  pp.  723-7  28 .  ~ 
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FIGURE  4  ENERGY  CURVES  FOR  NO* 
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FIGURE  6  TWO  INTERACTING  NO+  MOLECULES 
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FIGURE  6  ENERGY  CURVES  FOR  BROADSIDE  COLLISION 
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FIGURE  7  ENERGY  CURVES  FOR  HEAD-ON  COLLISION 
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The  ground  and  lowest  excited  states  calculated  for  the  broad¬ 
side  collision  are  given  in  Figure  6.  These  curves  were  obtained 
by  varying  the  N-0  bond  distances  together  and  by  the  same  amount 
in  both  molecules,  and  represent  the  effective  energy  per  molecule 
determined  by  dividing  the  energy  of  (N0+)2  by  2.  The  ground  state 
curve  <?Q°(2G)  is  somewhat  broader  than  the  1G  curve  for  N0+  in 
Figure  4.  Inasmuch  as  the  curves  approximate  harmonic  oscillator 
potentials,  the  frequency  relationship 


v 


2fB° 


SR' 


(46) 


may  be  used  to  show  that  the  low-lying  vibrational  energy  levels 
for  curve  2G  are  closer  together  than  those  for  curve  1G,  since 

32£?°(2G)  32£_B0(1G) 

- u  -  <  - V  -  #  (47) 

3R'  3R' 

Consequently,  compression  shifts  infrared  (IR)  spectra  to  lower 
frequencies.  If  the  NO+  molecules  were  pushed  closer  together  to 
a  separation  of  R  -  1.1&,  bonding  would  occur  between  nitrogens  to 
form  N2  and  between  oxygens  to  form  02+2. 

Compression  also  shifts  ultraviolet  (UV)  and  visible  spectra 
to  lower  frequencies.  This  is  seen  in  Figure  6  by  comparing  the 
curves  2E  and  2E '  at  R  =  4&.  The  allowed  vertical  transition  from 
2G  to  2E '  is  shifted  by  3.5  ev. 


Similar  shifts  of  UV  and  IR  spectra  are  found  for  the  head-on 
collision.  Figure  7.  These  curves  were  obtained  by  varying  the  N-0 
bond  distance  in  one  N0+  molecule,  and  represent  the  effective 
energy  per  molecule  determined  as  in  Figure  6.  The  curves  describe 
a  possible  situation  where  only  one  N-0  bond  is  breaking  during 
detonation.  The  ground  state  curves  4G  and  5G  indicate  substantial 
broadening  at  R  =  1.5&  and  spontaneous  decomposition  of  N0+  at 
R  =  1.15&.  Calculations  for  the  variation  of  two  N0+  bonds  yield 
similar  results.  Curve  3E  at  R  =  represents  the  same  state  as 
curve  IE '  in  Figure  4  for  the  isolated  molecule,  but  substantially 
lower-lying.  The  dashed  curve  3E'  at  R  =  2%.  represents  a  possible 
position  of  curve  3E  when  five  NO+  molecules  are  interacting  in  the 
linear  configuration  indicated  in  Figure  2.  This  curve  was  not 
calculated.  Curve  3E"  at  R  =  2&  is  a  dissociative  excited  state. 

It  intersects  the  3G  curve  at  2.5  bohrs. 


Figure  8  describes  a  possible  result  of  the  application  of 
degenerate  perturbation  theory,  using  the  perturbations  of  Equation 
(18)  and  the  matrix  element  in  Equation  (42) .  The  dashed  portions 
of  the  curves  were  not  calculated.  The  solid  portions  are  the 
curves  3G  and  3E"  of  Figure  7.  Note  that  a  new  activation  barrier 
Eact  woul^  be  created.  This  barrier  would  be  the  one  approximated 
by  Equation  (43) .  Calculation  would  show  the  size  of  the  gap  at 
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FIGURE  8  ENERGY  CURVES  RESULTING  FROM  COUPLING  OF  ELECTRONIC  AND  NUCLEAR  MOTIONS 
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the  top  of  the  barrier  to  be  2|£12|.  As  | ^12 1  increases,  so  does 
the  probability  P',  Equation  (39),  that  the  reaction  takes  place  on 
the  lower  curve.  Thus,  the  N-0  bond  io  N0+  will  break  more  easily 
when  (£12!  is  sufficiently  large  and  e|{1|  is  sufficiently  small. 

The  quantity  S12  *-n  Equation  (42)  will  be  calculated  and  the  results 
presented  in  a  future  report. 


The  curves  3G  and  3E"  are  electronic  states  of  different  symmetry, 
and  are  designated  as  i£+  and  ^tt,  respectively.  The  designation 
describes  a  spin  singlet  state,  S^q  =  0,  which  has  a  projection  of 
zero  electronic  orbital  angular  momentum  along  the  bond  axis,  Lz4'o  =  0. 
The  designation  describes  a  triplet  state,  S ^ H' ^  which 

satisfies  =  ±^1.  For  these  states,  <y0 1  He  |  «t'1>  =  0,  so  that  the 

conditions  in  Equation  (29)  are  satisfied  at  the  intersection. 


Figure  9  contains  the  repulsion  energies  for  the  (NO+) 2  system 
in  the  head-on  and  broadside  collisions.  The  estimated  pressures 
along  the  0°K  isotherm  are  indicated  at  various  intermolecular 
separations.  The  N-0  bonds  are  fixed  at  R^in  =  1.062&  for  NO  . 
These  pressures  were  calculated  from 


P  = 
o 


3R 


3£0 

9R 


cM 


(48) 


where  R^  is  the  distance  between  the  nuclear  centers-of-mass  of  the 
NO+  molecules.  Sq  is  the  effective  energy  per  molecule,  and  v  is 
assumed  to  be  the  average  volume  per  molecule,  v  =  R^,  which 

pertains  to  a  simple  cubic  lattice.  The  derivatives  3(?o/9R  were 
evaluated  numerically,  using  a  three-point,  divided-difference 
formula.  The  two  curves  do  not  represent  the  same  compression  since 
v  refers  to  the  separations  of  the  nuclear  centers-of-mass.  For 
example,  R^m  =  R  =  2%  for  the  broadside  collision,  and  RCM  =  3.062^ 
at  R  =  2%  for  the  head-on  collision. 


It  was  mentioned  in  the  previous  paragraph  that  the  derivative 
in  Equation  (48)  was  calculated  numerically.  This  procedure  is 

considerably  more  acurate  than  the  evaluation  of  from  the 

virial  theorem,  Equation  (12).  Although  Ke  is  calculated  in  the 
LCAO-MO-SCF  method,  it  is  not  determined  accurately  enough  relative 
to  So  and  9<?q/9R,  since  Equation  (12)  results  from  a  variational 
calculation  on  all  of  the  parameters  contained  in  the  wave  function 
To*  As  discussed  in  Appendix  A,  the  LCAO-MO-SCF  method  requires 
variation  of  only  linear  expansion  coefficients  in  the  MO's. 
Consequently,  Equation  (12)  applies  to  an  exact  Hartree-Fock 
calculation,  and  only  approximately  to  the  calculations  performed 
here. 
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The  pressures  shown  in  Figure  9  are  upper  limits  to  actual 
values  of  Pg  for  (N0+)2.  This  is  due  to  the  nature  of  the  methods 
used  here  and  of  the  exact  Hartree-Fock  methods  (HFM) .  Only  at 
very  high  compressions,  where  the  nuclear  repulsion  energy  dominates 
£q  ,  are  -£he  values  for  Po  accurately  determined.  At  low  compressions 
(R  >  2.5A),  configuration  interaction  (Cl)  methods  which  introduce 
electron  correlation  are  essential.  In  fact,  since  Hartree-Fock 
theory  does  not  introduce  correlation  between  electrons  of  the 
opposite  spin,  it  cannot  produce  a  minimum  in  the  interaction  energy 
corresponding  to  the  van  der  Waals  forces  between  molecules.  The 
most  important  term  due  to  correlation  is  the  dispersion  energy  Ed, 
defined  by23 


Ed  -  Eint  -  £hf  [<N0+>  2]  +  '  <« 

where  Ein{-  is  the  exact  interaction  energy  obtained  from  the  eigen¬ 
value  E  in  Equation  (24)  ,  £hf  [(N0+)  2]  is  Hartree-Fock  energy 

of  the  (N0+) 2  system,  and  £Hp(NO+)  is  the  Hartree-Fock  energy  of  the 
isolated  N0+  molecule.  In  the  LCAO-MO-SCF  method,  it  is  assumed 
that  the  errors  generated  by  not  performing  exact  Hartree-Fock 
calculations  are  similar  in  N0+  and  (N0+) 2  so  that  substantial 
cancellation  occurs  in  Equation  (49)  .  Then,  the  difference  in  the 
Hartree-Fock  energies  yields  the  classical  Coulomb  interaction 
energy,  the  first-order  exchange,  the  induction  energies,  some  of 
the  second-order  exchange  including  the  charge  transfer  part,  and 
some  higher-order  terms.  Reference  23  may  be  consulted  for  a 
discussion  of  these  contributions.  Since  Erl  is  a  negative  quantity 
(for  molecular  attraction)  which  varies  approximately  as  -R-6, 
the  derivative  of  Equation  (49)  shows  that  the  Pg  values  in  Figure  9 
are  too  high  by  ( 1/3R^M) (9E^/3R) . 

The  neglect  of  electron  correlation  also  leads  to  errors  in 
the  relative  placements  of  the  ground  and  excited  state  curves  in 
Figures  4-8.  The  results  of  the  CNDO/S-CI  calculations  show  that 
the  introduction  of  electron  correlation  can  reduce  the  separation 
of  the  curves  by  as  much  as  45  kcal/mole  for  (N0+) 2  or  23  kcal/mole 

for  N0+.  Consequently,  the  estimated  barrier  and  Eg^t  in  Fi9ure  8 
will  be  smaller. 
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SUMMARY  AND  FUTURE  WORK  TO  BE  PERFORMED 


The  theoretical  basis  for  an  analysis  of  ground  and  excited 
state  potential  energy  surfaces  has  been  presented.  The  behavior 
of  the  ground  and  excited  state  surfaces  of  molecular  systems  under 
compression  was  examined  within  the  framework  of  the  Born-Oppenheimer 
(BO)  approximation,  and  from  the  standpoint  of  intersections  of  the 
surfaces  which  result  in  a  breakdown  of  the  BO  approximation. 

An  application  of  the  theory  was  made  to  the  (N0+)0  system. 
Excitation  energies  Ae  for  the  first  excited  state  of  linear 
(N0+)n  decreased  from  Ae  =173  kcal/mole  for  n  =  1  to  125  kcal/mole 
for  n  =  2*,  and  to  Ae  =  100  kcal/mole  for  n  =  3.  The  molecules  were 
separated  by  3&  between  their  center s-of-mass.  For  3  <  n  <  6,  Ae 
increased  and  approached  a  plateau  at  -120  kcal/mole.  This  plateau 
is  close  to  the  value  of  a£  for  (NO+>2  which  suggests  that  the 
detailed  study  of  excitation  energies  and  activation  energies  can 
be  carried  out  on  just  two  or  three  interacting  molecules.  Minima 
were  reached  at  n  =  3  for  other  orientations  of  the  NO+  molecules. 
Recent  results  for  two  interacting  nitromethane  molecules  indicate 
that  a  similar  precipitous  drop  in  Ae  occurs.  Calculations  on 
nitromethane  are  continuing. 

It  was  seen  that  the  ground  state  curves  for  the  N-0  bond  stretch 
in  (NO+) 2  describe  a  bound  system  for  intermolecular  distances  greater 
than  about  1.15a.  At  distances  less  than  2&,  the  curves  are  broadened 
considerably,  indicating  a  weakened  bond.  At  the  distance  of  1.15&, 
the  curve  collapses,  suggesting  spontaneous  decomposition  of  NO+  at 
high  compression.  A  broadening  of  a  ground  state  curve  indicates 
that  the  vibrational  spectra  are  shifted  to  lower  frequencies.  The 
UV  spectra  are  also  shifted  to  lower  frequencies. 

Intersections  of  the  ground  and  excited  state  curves  are  found 
for  (NO+) 2  when  the  centers-of-mass  of  the  molecules  are  separated 
by  3&  or  less.  These  intersections  generate  new  barriers  to  bond 
scission  which  are  substantially  less  than  the  barrier  height  of 
-250  kcal/mole  in  the  gas  phase. 

These  calculations  on  (N0+)n  constitute  the  first  stage  in  a 
complete  series  of  calculations  which  will  also  involve  the  use  of 
perturbation  theory  to  lift  the  degeneracies  (intersections).  Some 
additional  mathematical  details  are  to  be  worked  out  in  the 
perturbation  analysis.  For  example,  the  wave  equations  for  the 
n-molecule  system  must  be  derived.  A  new  stage  of  calculations 
involving  these  wave  equations  will  complete  the  description  of  the 
kinetics  of  the  bond  scission  process,  and  will  be  presented  in  a 
future  report.  The  calculations  in  the  present  report  will  be 
refined  further  by  complete  Hartree-Fock  and  Cl  calculations  on 
the  excited  states.  The  sequence  of  calculations  on  (NO+)n  will 
serve  as  a  guide  for  those  to  be  done  on  (NO)n  and  nitromethane. 
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Appendix  A 


DESCRIPTION  OF  COMPUTATIONAL  METHODS 


The  treatment  of  the  electron  many-body  problem  in 
the  solution  of  Equation  (24)  is  usually  made  by  the  method 
of  configuration  interaction  (Cl).  Here,  the  eigenfunction 
is  expressed  as  a  linear  combination  of  Slater  determinants  Y j , 


00 


(A-l) 


Each  determinant  Yj  of  dimension  N,  the  number  of  electrons, 
is  constructed  from  N  spin-orbitals  selected  from  a  complete 
set  of  an  infinite  number  of  orthonormal  spin-orbitals, 

{uk}.  Each  spin-orbital  uk  =  pka  is  a  function  of  the 
coordinates  of  a  single  electron  so  that  Yj  can  be  expressed 
as 


Yj  =  a|Ui(1)u2(2) 


(A-2) 


where  A  is  an  antisymmetrizer  operator.  The  complete  set 
of  spin-orbitals  is  generated  from  a  Hartree-Fock  calculation 
which  minimizes  the  energy  given  by 

£0  "  <VHelV  '  (A“3) 

subject  to  constraints  on  the  orthonormality  of  the  spatial 
parts  pk  of  the  spin-orbitals.  The  solution  of  the  simultaneous 
equations 

X)  {tjIH.It,  -  EkSij}ckj  >  0  (A-4) 

yields  the  manifold  of  electronic  ground  and  excited  states. 
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The  solution  of  Equation  (A-3)  is  made  tractible  by  the 
careful  selection  of  combinations  of  spin-orbitals  which 
provide  the  most  rapid  convergence  in  Equation  (A-l) .  The 
success  of  this  procedure  depends  in  large  part  on  the  solution 
of  Equation  (A-3) ,  which  is  made  tractible  by  a  method  devised 
by  Roothaan  in  Reference  38.  This  method  is  called  the 
linear  combination  of  atomic  orbitals  (LCAO)  molecular  orbital 
(MO)  self-consistent  field  (SCF)  method,  the  LCAO-MO-SCF  method. 
Here,  each  MO  Uu  is  expressed  in  terms  of  the  finite  set 
of  AO's  { g i } r 


u 


k 


ckigi 


(A-5) 


A  careful  selection  of  a  basis  set  of  AO's  { g £ }  often  permits 
the  calculation  of  the  ground  state  energy  Cq  which  is  95-99% 
of  the  experimental  energy  E0  at  certain  nuclear  configurations, 
including  the  equilibrium  positions  of  the  nuclei.  As  a 
consequence,  Equation  (A-l)  reduces  to  4>q  =  Tq  for  a  closed- 
shell  molecular  ground  state.  In  the  case  of  an  excited  state 
corresponding  to  an  open-shell  electronic  configuration,  is 
often  expressed  in  terms  of  two  or  more  determinants.  The 
quantity  Er  -  <?K  is  often  referred  to  as  the  correlation  energy. 
This  energy  reflects  the  detailed  interactions  between  electrons 
which  were  ignored  in  the  solution  of  Equation  (A-3) . 

The  coefficients  c^j  in  Equation  (A-5)  are  obtained  by 
solving  the  pseudo-eigenvalue  SCF  equations 


A  A  AAA 

F  c  =  e  S  c  (A-6) 


in  conjunction  with  the  minimization  of  the  energy 

A  j  A  A  l  A  A  i  A 

£k  *  HTDt  +  IsdJp  -  JjdJq  .  (A-7) 

The  solution  of  these  equations  constitutes  the  variational^ 
method  of  the  SCF  approximation.  In  these  equations  F,  H,  P,  and 
Q  are  matrices  with  matrix  elements  containing  kinetic  energy, 
electron-nuclear  attraction,  and  Coulomb  and  exchange  electron- 
repulsion  integrals  over  the  AO's.  DT  and  Dq  are  density 
matrices  containing  the  coefficients  c^.  S  is  an  overlap 
matrix,  and  8  is  a  diagonal  matrix  which  contains  the  orbital 
energies.  Equations  (A-6)  and  (A-7)  frequently  provide  accurate 
descriptions  of  quantum  mechanical  systems,  and  can  usually 
be  solved  at  relatively  low  cost.  Great  strides  have  been 
made  toward  large  reductions  in  the  cost  of  these  calculations 
on  large  molecular  systems. 
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Equation  (A-7)  describes  the  averaged  interactions 
of  electrons.  As  a  consequence,  it  may  yield  surfaces  which 
have  incorrect  shapes.  This  problem  will  arise  when  the 
correlation  energy  varies  greatly  with  nuclear  configuration. 
Therefore,  a  compromise  must  be  struck  between  the  solution 
of  the  higher  cost,  accurate  Cl  methods  and  the  lower  cost, 
less  accurate  SCF  methods.  The  trick  is  to  reduce  the  number 
of  basis  functions  in  Equation  (A-l)  to  a  point  where  Cl 
methods  can  be  used  in  large  molecular  systems.  Then,  a 
variational  calculation  on  the  coefficients  DKj  in  a  limited 
basis  set  yields  a  finite  set  of  equations.  Equation  (A-4) , 
which  can  be  solved.  Perhaps  the  most  promising  approach 
introduced  so  far  is  the  Cl-related  method  called  the  generalized 
valence  bond  (GVB)  method. 18 

All  of  the  methods  described  above  are  ab  ini tio  methods. 
They  require  the  processing  of  M^/8  two-electron  integrals 
where  M  is  the  number  of  AO's.  A  great  deal  of  work  has 
gone  into  devising  semiempi rical  techniques  which  can  be 
used  for  large  molecules  at  low  cost.  The  one  which  is 
most  widely  used  is  the  complete  neglect  of  differential 
overlap  (CNDO)  method.  Here,  a  large  number  of  electron 
repulsion  integrals  are  eliminated  by  virtue  of  their  small 
values  as  dictated  by  the  magnitude  ofthe  differential  overlap 
of  AO's  occurring  in  them.  Also,  only  the  valence  orbitals 
are  treated  explicitly  in  the  SCF  procedure  while  the  core 
shell  contributions  are  fixed  and  determined  by  parameters 
defined  in  terms  of  experimental  ionization  potentials. 

Parameter izations  in  the  method  have  been  refined  to  the 
point  where  reliable  results  can  be  obtained  for  a  wide 
variety  of  molecular  properties.  A  method  which  yields 
reliable  electron  excitation  energies  for  many  systems  by 
combining  the  best  features  of  the  CNDO  and  Cl  techniques 
is  called  the  CNDO/S-CI  method. 39  it  employs  equations 
similar  to  Equations  (A-l)  -  (A-7) .  This  technique  is  of 
particular  value  for  the  present  calculations. 
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